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Abstract 



We have investigated the quantum-to-thermal crossover temperature T* 

for cavitation in liquid helium mixtures up to 5% 3 He concentrations. With 

respect to the pure 4 He case, T* is sizeably reduced, to a value below 50 mK 

for 3 He concentrations above 2%. As in pure 4 He, the homogeneous cavitation 

pressure is systematically found close to the spinodal pressure. 
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Quantum cavitation in superfluid liquid 4 He has been experimentally observed by Balibar 
and coworkers [jXJ] . They have used an ultrasound technique by means of which, a large 
pressure oscillation is created in a small bulk region of the experimental sample. This 
prevents cavitation on the walls of the experimental cell. Moreover, as superfluid 4 He can 
be made perfectly pure, heterogeneous cavitation can be avoided. Due to the nature of 
these experiments, it is very difficult to determine the pressure (P) and temperature (T) at 
the cavitation site. Consequently, a quantitative analysis of the experimental results relies 
on theoretical estimates of the equation of state and cavitation barriers in the spinodal 
region |||| . Depending on these theoretical inputs and on experimental conditions, such as 
volume and observation time, it is inferred that quantum cavitation likely takes over thermal 
cavitation at a temperature in the 200-240 mK range. The agreement between experiment 
lU and theory can be considered as satisfactory given the inherent limitations of both. 

In this work we address the problem of thermally assisted quantum cavitation in low 3 He- 
concentration, liquid helium mixtures (up to 5%). The relevance of this problem is twofold. 
First, the planned experiments on quantum cavitation in these mixtures ||. Second, the 
complexity introduced in the theoretical description of the still superfluid liquid 4 He, when 
a small but sizeable 3 He amount is present. To our knowledge, no investigation of this kind, 
either experimental or theoretical, exists in the literature. Some effort has been concentrated 
in the study of supersaturated 3 He- 4 He mixtures at positive pressures and low temperatures 
(see for example Refs. 0-0), or near the tricritical point (see Ref. || and refs. therein). 

The method we use here is conceptually simple and technically workable. It is based, 
on the one hand, in using a density functional [|1(J that reproduces the thermodynamical 



properties of the 3 He- 4 He liquid mixture at zero temperature. This functional has been 
slightly modified to reproduce the surface properties of the 3 He- 4 He interface, and those 
of the mixture free-surface |TI],fnj. All these properties are relevant for a quantitative de- 
scription of the cavitation process. On the other hand, we make use of a functional-integral 
approach especially well suited to find the crossover temperature T* with the only further 
approximation of imposing irrotational motion for the growing of the critical bubbles. 



Let us first generalize the functional-integral approach as used in 0, to the mixture 
case (see also |13|] for another example of how it applies). We recall that above a crossover 
temperature T*, the cavitation rate, i.e., the number of bubbles formed per unit time and 
volume, is given by 

J T = J 0Te - An — / T , (1) 

where AQ max is the barrier height for thermal activation. Below T*, assuming AQ max » 
T, one can write 

J Q = J 0Q e- s ™, (2) 

being SVmn the minimum of the imaginary-time action 

S(T) =j&rj&r C. (3) 

£ is the imaginary-time classical Lagrangian density of the system. For a given T, the 
time-integration is extended over a "period" r in the inverted barrier potential well with 
t = h/T. The crossover temperature T* is determined from the angular frequency u p of the 
small amplitude oscillations around the minimum of the potential well (see |TJj and Refs. 



therein). The prefactors Jot and Jog depend on the dynamics of the cavitation process. 
In actual calculations, they are often estimated as the number of cavitation sites per unit 
volume times an attempting frequency. 

The Lagrangian density can be easily obtained from the zero-temperature density func- 
tional of [pT|,r2[ , whose use is justified in view of the low temperatures involved (below 



200 mK). The critical cavity density profiles p°(r), p°(r) are obtained solving the coupled 
Euler-Lagrange equations 

= ,g = 3,4, (4) 

where w(p3, p±) is the grand potential density and p q are the particle densities of each helium 
isotope. Afl max is given by 



Aa 



dr 



^(P^Pa) - ^{Pm3,PmA) 



(5) 



where p mq are the corresponding densities of the metastable homogeneous liquid (see [flT|J12 
for details). Assuming that only spherical bubbles are developing, the collective velocities 
t) of both helium fluids are irrotational and one can then define for each isotope a 
velocity potencial field s q (r, t) such that u q (f,t) = V s 9 (r, t). It follows that 



£ = m qPq s q ~ ^G°3, Pa, s 3 , Sa) , 



(6) 



where Ti.(p q , s q ) is the imaginary-time Hamiltonian density 



^(P3, PA, S 3 , Sa) = ~Y1 m qPqUq 2 ~ [ UJ {P3, Pa) ~ ^(Pm3, Pwa)\ ■ 

1 q 

Hamilton's equations yield the following four equations : 

5H 



m q p q 



Ss, 



-m q V(p q u q 



(7) 



(8) 



m q s q 



5p q 



Eqs. (H) are the continuity equations. Taking the gradient of Eqs. (||) one gets 



dug 
dt 



To determine T* one has to find the small amplitude, periodic solutions of Eqs. 
(0) linearized around p° and p° f3|,|l3|]. Defining the transition densities pl(r) as 



(9) 



(10) 
and 



Pq(r, t) = p° q {r) + p\(r) cos{u p t) 



(11) 



where p q {r) is much smaller than p q (r), and keeping only first order terms in u q (r,t) and in 



pUr), one gets: 



u 2 p p\{r) 



P° q (r)V £ 



5 2 u 



p l Ar) 



9 = 3,4. 



(12) 



In Eq. 



5 2 u . A 



Sp q Sp q 



- • p l q i{r) means that 5uj/5p q has to be linearized, keeping only terms in 



p\ and p\, and their derivatives. 



Eq. (|12|) is a fourth-order linear differential, eigenvalue equation for the "vector" 
(pl(r), p\ (r)), whose right-hand side term is straigthforward but very cumbersome to obtain. 
We have done it using the Mathematica software package ITS] . As in |J, a careful analysis 
shows that the physical solutions of Eqs. ( |T2"D have to fulfill (Pq)'(O) = (Pg)"'(0) = 0, and 
have to fall exponentially to zero at large distances. From the linearized continuity equation 
p l q {r) oc — V(pqMg), it is obvious that the integral of p^ir) over the whole space is zero. 

We have solved the eigenvalue Eq. (|l^) as in |§. For a given pressure and 3 He- 
concentration x = Pm3/(p m 3 + Pm*), om y a positive eigenvalue is found, out of which 
we get T* = hujp/27i. 

Fig. 1 shows T* (niK) as a function of P(bar) for x = 0.1, 1, 2, 3, 4 and 5 %. Compared 
to the pure 4 He case [[J, T*(P) has now a more complex structure. It is worth to note that 
the maximum of the T*(P) curve has decreased from ~240 mK for pure 4 He || down to ~ 
140 mK for 3 He-concentrations as small as 1%. 

Fig. 2 shows two different bubble configurations for x= 1%. Configuration (a) corre- 
sponds to P= —8.23 bar and T* = 62.7 mK, and configuration (b) to P= —5.17 bar and 
T* = 115.7 mK. The solid lines represent the 3 He and 4 He critical bubble particle densities 
in A , and the dashed (dash-dotted) lines represent pl(r) (p\(r)) in arbitrary units. Near 
the spinodal region, the "bubble" configuration is filled with 3 He: the surface tension that 
matters for bubble formation is that of the 3 He- 4 He interface. Away from the spinodal region 
(configuration (b)), the critical bubble is a true bubble covered with 3 He (Andreev states): 
the surface tension that matters now is that of the 3 He- 4 He liquid free-surface, which is 
about ten times larger than the previous one. The different surface tensions involved in 
these processes, together with the existence of a 3 He- 4 He saturation curve at negative pres- 



sures down to x ~ 2.4% |TT| are the cause of the structures displayed in Fig. 1. We will not 
give here any further detail since, as in the pure 4 He case, only the part of the T*(P) curve 
near the spinodal region is relevant for the cavitation problem. 

It is interesting to see that the transition densities p 1 evolve from those corresponding to 
"volume oscillations" (Fig. 2a) to "surface oscillations" for 4 He, and a mixed surface- volume 



type for 3 He (Fig. 2b), to eventually become pure "surface oscillations" for both isotopes 
when moving from the spinodal towards the saturation line. This fact has been already 
found in pure 4 He ||. 

To establish which T*(P) in Fig. 1 corresponds to that of the experimental conditions, 
we have to determine the homogeneous cavitation pressure P^, which is the one the system 
can sustain before bubbles nucleate at an appreciable rate. P^ can be obtained solving the 
equation 

l = (Vt) e J, (13) 

taking J = Jy and Jot — {ksT) j '(hVo) if T>T*, where ksT/h is a thermal attempting 
frequency, and Vq = AtcRI/3 is the inverse of the number of possible nucleation centers per 
unit volume, with R c =10 A as a typical radius of the critical bubble. For T < T*, one 
takes J = Jq and, lacking of a better choice, Joq = Jot(T=T*). It yields Joq of the order 
of 10 5 A 3 s _1 . Another possible guess for the quantum attempting frequency is to equal it 
to v s /R c , where v s is the sound velocity, which is of the order of 10 12 A s _1 in the region 
of interest. This choice differs from the previous one in two orders of magnitude. However, 
this difference in the prefactor value does not alter in practice the results. Indeed, we have 
chosen two extreme experimental values for the factor (Vt) e (experimental volume x time) 
which differ ten orders of magnitude without much influencing the obtained T*-value. 

P^(T*) is shown as circles (squares) on the curves in Fig 1. The circles correspond 
to (Vt) e = 10 4 A s, and the squares to 10 14 A s. Compared to the pure 4 He case [|[, and 
depending on the (Vt) e -value, for x = 1% T* has been reduced a factor of 4 or 5, respectively. 

2 3 

Fig. 3 shows Ph as a function of T for the above mentioned ^-values and (Vt) e = 10 14 A 
s. Thermal and quantum regimes are displayed. The dashed line is the extrapolation of the 
thermal regime to temperatures close to T=0. Finally, Fig. 4 displayes the T=0, P sp (x) 
spinodal line together with the P/ l (x,T*) curves for the indicated (Vt) e values. Notice that 
the smallest x-value displayed in this figure is 0.1%. For 3 He concentrations closer to zero, 
it is very unlike that 3 He has time enough to difuse and develop the critical configurations 
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that constitute the starting point of the present calculations. Also shown in that figure (dots 
on the y-axis) are the pure 4 He values || . 

In conclusion, we have thoroughly investigated thermally-assisted quantum cavitation 
in 3 He- 4 He liquid mixtures at small 3 He concentrations. Our quantitative predictions rely 
on a robust density functional that reproduces the relevant characteristics of the mixture 
and its interfaces, and on a sound formalism, the functional integral theory. Besides the 
approximations inherent to the method (use of a density functional, irrotational flows), 
no further approximation has been introduced to obtain the quantum-to-thermal crossover 
temperature. The present results might thus be a valuable guide to the planned experiments 
I- 
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FIGURES 

FIG. 1. T* as a function of P. The homogeneous cavitation pressure Pft(T*) is shown as circles 
(squares) for (Vt) e = 10 4 A 3 s (10 14 A 3 s), respectively. 

FIG. 2. Particle densities p\{r) and p%{r) of the critical bubbles (solid lines), and the p\{r) 
(dash-dotted lines) and p\{r) (dashed lines) transition densities for x= 1%, corresponding to: (a) 
P= -8.23 bar and T* = 62.7 mK. (b) P= -5.17 bar and T* = 115.7 mK. p\{r) are drawn in 
arbitrary units, p° q {r) in A -3 . 

FIG. 3. Homogeneous cavitation pressure P/j as a function of T for the same x-values as in 
Fig. 1. 

FIG. 4. The T=0 spinodal line P sp (x) (solid line), and the homogeneous cavitation pressure 
P/ l (x,T*) (dashed and dash-dotted curves) for the indicated (Vt) e values. The dots on the y-axis 
are the corresponding pure 4 He values. 
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